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Abstract 

The renormalization factors of the dimension-six effective operators for proton 
decay are evaluated at two-loop level in the supersymmetric grand unified theories. 
For this purpose, we use the previous results in which the quantum corrections to 
the effective Kahler potential are evaluated at two-loop level. Numerical values for 
the factors are presented in the case of the minimal supersymmetric SU(5) grand 
unified model. We also derive a simple formula for the one-loop renormalization 
factors for any higher-dimensional operators in the Kahler potential, assuming that 
they are induced by the gauge interactions. 



1 Introduction 



Discovery of the Higgs boson [H[2] may suggest the existence of supersymmetry (SUSY). 
The supersymmetric theories may accommodate the hierarchical structure with the great 
desert naturally. Searches for rare processes, such as proton decay, would be one of meth- 
ods to access the physics beyond the supersymmetric standard models (SUSY SMs). The 
processes are dictated with the effective higher- dimensional operators. When comparing 
the prediction with the observation precisely, we need to include the radiative corrections 
correctly. 

The realization of the gauge coupling unification strongly motivates us to study the 
supersymmetric grand unified theories (SUSY GUTs) [3H5] . In the theories, proton decay 
is induced by the exchanges of the colored Higgs multiplets and the X gauge bosons, which 
yield the baryon and lepton number non- conserving interactions. They are expressed in 
terms of the dimension-five and -six effective operators, respectively. It is found that the 
former interactions in general give rise to dominant channels for proton decay, such as 
p — > K + u [HE]- However, the current experimental limits on the channel, r(p — > K + v) > 
3.3 x 10 33 yrs [9], are so severe that the contribution of the dimension-five operators is 
required to be suppressed by a certain mechanism; otherwise the model is excluded just 
as the case of the minimal SUSY SU(5) GUT unless the SUSY particles in the SUSY 
SM are much heavier than the weak scale [10, llj. A variety of such mechanisms have 
been proposed. For example, the Peccei-Quinn symmetry [T2] would be exploited for 
the purpose. The R symmetry also plays a role in suppressing the dimension-five proton 
decay in the models with extra dimensions [13]. With such a suppression mechanism 
imposed, the dimension-six operators in turn become dominant. In this case, the main 
decay mode is the p — > 7r°e + channel; the present experimental limit on its lifetime is given 
by r(p ->■ 7r°e + ) > 1.29 x 10 34 yrs [H]. Since in the SUSY GUTs the GUT scale M GUT is 
relatively high, i.e., M GVT ~ 2 x 10 16 GeV, the predicted proton lifetime usually evades 
the experimental limit. However, the consequence might be altered if there exist extra 
particles in the intermediate scale. With such particles belonging to a representation of 
the grand unified group, the gauge coupling unification is still achieved, while its value 
at the unified scale turns out to be enhanced. Then, the proton lifetime is considerably 
reduced due to the large gauge coupling [To] . 

In order to study such possibilities based on the proton decay experiments, it is impor- 
tant to make a precise prediction for the decay rate. To that end, we need to determine 
the effects of the dimension-six operators, which are generated at the GUT scale, on the 
low-energy physics by using the renormalization group equations (RGEs). Indeed, there 
have been several literature in which the renormalization factors for the effective operators 
are evaluated. In Ref. [IB], the long-distance QCD corrections are computed at two-loop 
level. For the short-distance factors, on the other hand, only the one-loop calculation is 
carried out in Ref. p2J in the SUSY SM. 

In this Letter, therefore, we evaluate the renormalization factors of the dimension-six 
operators at two-loop level in the presence of the supersymmetry. In the calculation, 
we use the results for the two-loop corrections to the effective Kahler potential given in 
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Ref. [18J. Since in the SUSY GUTs, the most of the intermediate energy scales are su- 
persymmetric, the short-distance renormalization factors are well approximated by those 
evaluated in purely SUSY theory. Thus, combined with the long-distance effects given in 
Ref. [16J, our results offer a tool for making a prediction of the proton decay rate with 
accuracy of two-loop level. 

We also derive a simple formula for the one-loop level renormalization factors of any 
higher-dimensional operators in the Kahler potential, when including only the gauge inter- 
action contributions. It is applicable to other observables, such as the neutron-antineutron 
oscillation [19]. 

This Letter is organized as follows: in Sec. [21 we first write down the dimension-six 
effective operators in terms of the superfield notation. Our notations and conventions 
are also shown in the section. Then, in the subsequent section, we describe a way of 
calculating the renormalization factors of the operators by using the effective Kahler 
potential, and present the results for the computation. In Sec. IU the comparison of the 
one- and two-loop renormalization factors is discussed in the minimal SUSY SU(5) GUT. 
Section [5] is devoted to conclusion and discussion. 



2 Dimension-six effective operators 

To begin with, we write the dimension-six effective operators for proton decay in a SUSY 
and gauge invariant manner with superspace notation: 

000 = J tfetffc^e^efrKie-^rfyQfal , (1) 

where all the chiral superfields correspond to left-haded fermions, and V\ and V3 are 
the U(l)y and SU(3)c vector superfields with the gauge coupling constants gy and g 3 , 
respectively. The subscripts i,j, are the SU(2)i indices, while a, b, c are the color indices. 
Furthermore, we omit the generation indices for simplicity. 

The relationship between bare and renormalized operators is written in the following 
form: 

Of = Z {1) {1 \ (1 = 1,2), (2) 

where the subscript B indicates the operator is bare. Then, the Wilson coefficients C^ 1 ' 
for the operators obey the differential equations, 

^C^^) = loW C^^), (3) 

with 7 C (j) the anomalous dimensions for the operators defined as 

700,=^ In Z«. (4) 
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The anomalous dimensions are obtained by analyzing the vertex functions (or the effec- 
tive action) in which the operators are inserted. Since we now deal with the dimension-six 
operators which contain four chiral or anti-chiral superfields, it is sufficient to consider the 
four-point vertex functions which include the corresponding external superfields. Their 
renormalization group equations (RGEs) are given as 



d a d 
dp dg 



5^ 7i + loc 



r 



. 



(5) 



Here, Tqw are the four-point vertex functions with an insertion of the operators O^ 1 '. The 
gauge coupling constants and their beta functions are denoted by g a and j3 a , respectively, 
and the sum over each gauge group is implicit. Further, ji shows the anomalous dimension 
of each superfield contained in the operators. From now on, we often omit the superscript 
(/) for brevity. 



3 Renormalization factors 

In this section, we present the formulae for the renormalization factors. They are derived 
from the effective Kahler potential given in Ref. [18]. In the calculation, the dimensional 
reduction scheme (DR) [20J is employed for the regularization. We first obtain the one- 
loop results and confirm the results in Ref. p2] in the former subsection. Then, in the 
latter subsection, we evaluate the two-loop contribution. 

3.1 One- loop 

Let us first evaluate the vertex functions at one-loop level. For this purpose, we use the 
results in Ref. [18], where the effective Kahler potential for generic four dimensional N — 1 
SUSY theories is computed up to two-loop level. According to the results, the one-loop 
correction to the Kahler potential is given as 

1 j^f2 



where /1 2 = Ane 7 /i 2 defines the MS renormalization scale, and the mass matrix Mq,s is 
defined by 

( 71,, 

L C( a ), 

with cf) the background for the chiral superfield $ and G a b the Kahler metric 



(Ml M ) AB = 2glUTt ) )\G\(T£ ) y i 4> d , (7) 



In Eq. (jSJ), Tr denotes the trace over the adjoint representation of a gauge group whose 
coupling constant is g a and generators are given by t\ . Moreover, in the following 
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calculation, we only take the gauge interactions into account, i.e., we neglect the super- 
potentialQ 

In order to obtain the renormalization factors for the higher- dimensional effective 
operators, we consider the Kahler potential 

K = a a + CO + CO ] , (9) 

with C the Wilson coefficient of the operator O. In this case, the Kahler metric reads 

G\ = 5 a b + CO\ + CO\ , (10) 

with O a b = d 2 (D/d<f) a d<f) b . By substituting the above equations to Eq. ([6]), we have 

AK ^ = ~ E + ^ 2 )[C a {a)M a + {C^T { A a) ) a O%(T^U) b + h.c.}] , (11) 

where 6*0,(2) are the quadratic Casimir group theory invariants for the superfield 
defined in terms of the Lie algebra generators T4 by (Tj^T^ )g = C a (i)5 b . Further, we 
keep only the terms up to the first order with respect to the Wilson coefficient, C, and do 
not show the terms including the logarithmic dependence on the background fields, which 
are not relevant to the present calculation. At the first order in the perturbation theory, 
the RGE (El) then leads to 

i a 

Here, the superscript (1) of the anomalous dimensions denotes that they are evaluated at 
one-loop level. In supersymmetric theories, 7^ is given as^| 

7 « = -2X^(0^. (13) 

a 

Now we evaluate the second term in Eq. (1T2|) . To that end, we analyze the structure of 
the term on a general basis in order to derive the formula for the one-loop renormalization 
factor of any operator. Consider the following operator which contains an arbitrary num- 
ber of both chiral and anti-chiral superfields and is singlet under a given global symmetry 
G as a whole: 

o = Kr im K-^ ■ ■ ■ ^ • <I>/; • • • 'I*" • (14) 



1 Experimental constraints on the effective operators in Eq. (TT]) are particularly severe when the 
external lines of the operators are of the first and/or second generations. In such a case, the size of the 
Yukawa couplings are negligible. 

2 The anomalous dimension of fields, ji, may be also derived at one- and two-loop levels from the 
effective Kahler potential derived in Ref. Q2] in the similar way. See the first term in Eq. ([TTj). 
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Here, the coefficients A" Jn and A^ 1 "- 1 ™ make the set of superfields G singlet. When G 
is localized (gauged), the operator invariant under both supersymmetry and the gauge 
symmetry is 

J d 2 9<?9 (Ajr*»¥ 4l . . . $J [e^ TA ] \{\ h h .., n ^ ...<!*), (15) 

where g and are the coupling constant and the gauge vector superfields of the gauge 
group G, respectively Moreover, T A are assumed to be the generators for an irreducible 
representation, which are relevant to the transformation properties of the composite chiral 
superfield ^ . . . $j n ; under the gauge transformation, ^ . . . $ Jn is transformed as 

(K... jn ® jl ■ ■ ■ (e i9AAT y b (*L.in® jl ■ ■ ■ ' l> " ) > (Hi) 

with A A any chiral superfields. Further, we write the generators for each chiral superfield 
$ as t A , i.e., $ J — > (e l9A 1 Y ,$ J . Then, since the coefficients A? and A* 1, Jm assemble 
the transformation properties of each chiral superfield into that of the composite operator 
XI 4 . . . it follows that 

(T A r^l., n = %..,„(* A )i + • • • + K.., n -^ A A > ( 17 ) 

and similarly for the anti-chiral superfields, 

A>"^(T A ) 6 a = (^)«, \f i2 - im + ■■■ + {t A )\ A* 1 -*™- 1 *™ . (18) 

These expressions imply that Jn and A^ 1 '"*" 1 are invariant tensors under G. By using 
the relations, we now evaluate the second term in Eq. f fT2|) . It goes as follows: 

(4>t A ) a O a b (t A <p) b = [(t A )\ 2 ~ im + ■■■ + {t A )\ a* 1 -*™- 1 ^] 

x [%,.,^ A )i + • • • + ^...^(^Ij^ • • • 4^ • • • 0^ 
= a;;--- ! 7- , ) , vv- , )'',a^.. ; ,o, . . . ^ . . . <p 

= C™ O, (19) 

where Cg° mp is defined by T A T A = C^ omp ll; it corresponds to the Casimir invariant for the 
composite chiral superfield A| x Jn $- 71 ...$ Jn . Substituting the expression into Eq. f lT2|) . 
we finally obtain a generic formula for the one-loop renormalization factors of arbitrary 
operators: 

,(l) _ 9o 



(20) 



with C*™ mp the Casimir invariants of the gauge group a for the chiral part of the operators. 

So far we have assumed that the set of chiral (anti-chiral) superfields forms an irre- 
ducible representation. When it is reducible, independent operators are formed. They 
are not mixed with each other at one-loop level if only gauge interactions are effective. 
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Now we apply the formula to the dimension- six effective operators for proton decay in 
Eq. (pp. We find C 3 comp = C 3 (D) = 4/3 in the case of SU(3) C and C 2 comp = in the case 
of SU(2)l for both and . Here, □ denotes the fundamental representation of the 
corresponding group, and we have used C^iD) = C^(D). Note that the latter equation 
for SU(2)l follows from the fact that the SU(2)l non-singlet superfields in the effective 
operators have the same chirality and form an SU(2)i singlet. For U(l)y contributions, 
on the other hand, we obtain different results for the operators and O^: Cy° mp = 
(Y Q + Y L f for £> (1) and C™ mp = (2Y Q ) 2 for £> (2) . As a result, by using these factors we 



obtain that 



where 



-,2 



a=Y,2,3 



[?&],= [■$>]. = -!. (22) 



3 ~~ hoi 2 )] 


3 


8 

~3 


2 = bo(2)] 


2 _ 


-3 


0(1) J Y 


11 




XD 1 _ _ 


23 
~9~ 





(24) 

These results are totally consistent with those in Ref. [T7] . 
3.2 Two-loop 

Next, we discuss the two-loop level contribution. Again, we use the results in Ref. [IS] . 
The radiative corrections to the Kahler potential at two-loop level are described by 

f\K — 1 R b d T ac (M 2 \ \^ f( a ) f( Q ) jBDEAFC (M 2 \ 
AA 2 - -K a C J h d\ M ) - JABCJDEF 1 \ M V{a)) 

a 

- J2( GT A^) b M T B a) G)a d HY d AB (M^ M* (a) ) , (25) 

a 

with /abq the structure constants of the gauge group a. The mass functions and the 
geometric factors appear in Eq. (1251) are displayed in Appendix. By using them, we 
readily obtain the two-loop corrections to the vertex functions. We found from explicit 
calculation that the two-loop correction is not given simply by the gauge transformation 
properties of the composite chiral superfield in the operator and anomalous dimension of 
the external fields, which is different from the one- loop ones. At present, however, since 
the explicit derivations are quite complicated, we simply give the final results and defer 
full details 
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The RGE in Eq. (j5J) at two-loop level is given as 

^ + E " E 7, (1, rg» - E % ,2, rS* + 7g'rg» + o^'rg' = o . (26) 

Here, the subscripts (0-2) indicate the quantities are evaluated at tree, one-loop, and 
two- loop level, respectively. One-loop anomalous dimensions 7}^ are shown in Eq. ( !T3|) . 
while the two-loop ones are given as [21] 

if ] = j^-Ty E + 2g}C p (i)] . (27) 

Here, b a are the one-loop beta function coefficients for gauge coupling constants, given 
as b a = E7.T a (z) — 3C a (G) with C a (G) and I a (i) the quadratic Casimir invariant for the 
adjoint representation of the group a and the Dynkin index of the chiral multiplet 
respectively. 

From the RGE in Eq. ( |26l) . we now obtain the two- loop anomalous dimensions for the 
effective operators. Again, we parametrize them as follows: 



,.(2) = 93 r (2) I 9$ r (2) 1 gy r (2) 1 

~ (167T 2 ) 2 ^ '° W -' 33 (167T 2 ) 2 L ^C' 1 '- 122 (1 6^-2^2 L W J 



7T" / 

, #2^3 r (2) 1 , 9y92 r (2) i ffyffl r (2) i ^ 

(167T 2 ) 2 ^ C)(/) -' 23 (167T 2 ) 2 ^ '0(- r )jy2 (167T 2 ) 2 ^ ^d (I) Jy3 ' V ' 

Then, we have 

bS\,] 33 =bS,] 33 = ^ + 86 3 , (29) 



YY 



bS)] 22 = [tS)] 22 =~ + 36 2 , (30) 
113 5 



[7gi)]yy=§ + 36K, (31) 



[tSo] 2 3 = 20 i (32) 



boci)] 23 - 12 
( 2 ) 1 

'0( 2 )J 2 3 

boa)]y 2 = 2 ' 

[7S)]y 2 = ^, ( 33 ) 

u( 2 ) 1 = — 

L'o(i)Jv3 g ' 

= (34 > 
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4 Results 



In this section, we give the numerical results of the renormalization factors in the minimal 
SUSY SU(5) GUT. The short-distance renormalization factors A§ are defined as the 
ratios of the coefficients C^ 1 ' for the effective operators at the SUSY scale Msusy to those 
at the GUT scale M Gl]T : 

where we assume Msusy = 1 TeV and Mqut = 1-5 x 10 16 GeV. The numerical results at 
one-loop level are given as 



(36) 



^(lloop) = 1.960 , 
A^lloop) = 2.059 , 

while at two-loop level, we have found 

4 1) (21oop) = 1.962, 
^ 2) (21oop) = 2.053 . 

Here, we calculate the one-loop (two-loop) short- distance factors with the one-loop (two- 
loop) renormalization equations for the gauge coupling constants in the SUSY SM [21] . 
There is a cancellation among the two-loop corrections since the signs of [7^)] 33 and 

[7^(2)] 33 are opposite to those of the other two-loop anomalous dimensions. Therefore, 
the numerical values at two-loop level hardly differ from the one-loop ones. 



5 Conclusion and discussion 

We have evaluated the short-distance renormalization factors for the dimension-six proton 
decay operators at two-loop level with the effective Kahler potential. The procedure 
described in this Letter is generic and applicable to any higher- dimensional operators. We 
get the results A^(21oop) = 1.962 and A^ 2) (21oop) = 2.053 in the minimal SUSY SU(5) 
GUT. We have found that the two-loop contributions hardly change the renormalization 
factors evaluated at one-loop level. 

Finally, we briefly comment on the extensions of the minimal SUSY GUT. The gauge 
coupling constants at the GUT scale increase if there exist extra particles in the inter- 
mediate scale. The two-loop effects may be more significant in such cases. In addition, 
let us note that our results are only for the SU(3)<7X SU(2) £ xU(l)y gauge interactions. 
If some new gauge interactions exist below the GUT scale, we also need to evaluate the 
contributions of the gauge interactions. Even for such theories, however, it is possible to 
execute the prescription describe above to estimate the renormalization factors by means 
of the effective Kahler potential. 
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In this Letter, we neglect the possible effects of the threshold corrections from particles 
whose masses are around the GUT scale. Although the effects are model-dependent, to 
complete the two-loop level calculation, we also need to evaluate such corrections. We 
will discuss the issue on another occasion |19|. 
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Appendix 

Here, we show the explicit form of the mass functions as well as the geometric factors 
given in Eq. fl25|) : 

JVM 2 ) = TT^O^ 2 ) T,( M " G ~ 1 ) a b(MlG- l r d , (38) 

a,P 



I ABCDEF m (a) ) = ~\^^(lnfi 2 )[4(M^ a) ) AB 5 CD 5 EF 
- 5 AB (M^ (a) In M^cdSef - S AB 5 CD (M^ (a) In M* (a) ) EF ] + cycl. , (39) 
where the -l-cycl. denotes the cyclic permutation of the labels AB, CD, EF, and 

^V/ S (M 2 ,M^ (Q) ) = --§-(ln^)x 



—l\c 

d 



IQtt 

£ 5 AB {2(M^G- 1 r b (G- l ) c d + 2(G- 1 y b (M*G- 1 

- (G-r^HM^G-^d- (MlHMl}G^r b (G-^ d } 

+ 2{G- l )\{G- l y d {M 2 v(a) ) AB + (G-'UG-^M^lnM^)^ 



(40) 



Here, we drop the terms independent of the scale fi or containing two logarithms. The 
latter terms give rise to the logarithmic terms after differentiation, which cancel other 
logarithmic terms in the RGEs. The mass parameters are defined as 

(M*)\ = 2gl(T A a UY^T A a ~>G) bl (41) 
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and 

{M 2 v(a) ) AB = i [(M 2 c(a) ) AB + [M 2 c{a) ) BA ] . (42) 

Further, G~ l is inverse of the Kahler metric G a b defined in Eq. (jSJ), and the curvature 
R a b c d is given by 

R\ c d = -J^T h G c d ~ (4rG C e) (G-'Yf (4r h G f d \ . (43) 



The third term in Eq. ( 1251) includes the shorthand notations, {GT A (p) b . c and (4>TbG) c 
which are defined as 

{GT A <j>)\ = G a c (T A ) c b + (4^G a b ) {T A <j>Y 



(T A ) a c G c b + (4>T A ) C ( 4rG\ ) = (4>T A G) b ' a . (11) 



M 

Here, the second line follows from the gauge invariance of the Kahler potential. 
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